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Abstract. We compute the anomalous dimensions of the transversity operator at three loops 
in the MS scheme for fixed moment n where n < 8. The results for the RI' renormalization 
scheme are also provided for an arbitrary linear covariant gauge for n < 7. 
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Recently, the three loop anomalous dimension for the unpolarized twist-2 flavour non-singlet 
and singlet Wilson operators have been determined analytically for all values of the operator 
moment n, [U El El S]. This was a formidable undertaking spanning ten years and relied on 
cutting edge computer algebra and symbolic manipulation techniques implemented on high 
performance computers. With the obvious necessity of such results to ensure the complete two 
loop evolution of the Wilson coefficients and obtain more precise estimates of quantities such as 
as{Mz), there is also the need to extend such computations to problems involving spin. Indeed 
in this area, one quantity which will be of interest is that relating to transversity which was 
originally introduced in [H (H [7]. This corresponds to the probability of finding a quark in a 
transversely polarized nucleon polarized parallel to the nucleon versus that of the nucleon in 
the antiparallel polarization. Although experimentally it is harder to extract information on the 
transversity compared with usual deep inelastic scattering, one will still require the anomalous 
dimensions of the underlying operator to as high a loop order as is calculationally feasible for 
accurate renormalization group evolution. Currently the two loop anomalous dimensions are 
available for arbitrary moment, n, in the MS scheme, [H O [TOl ITTl I12j . However, given the 
symbolic manipulation machinery now available [H [2l [3l HI [13] , it is clearly only a matter of time 
before the full n-dependent three loop MS results are determined. Although the underlying 
transversity operator differs from that of the Wilson operators of [B El [3l H] , calculationally it is 
on a par with the non-singlet work of [14j. Prior to the full n-dependent results of [H El Eld] for 
the twist-2 Wilson operators, one approach was to carry out a fixed moment determination of the 
anomalous dimensions. Essentially the first even moments to n = 16 (apart from n = 14) were 
determined, [HI |T5l [16] . Although eventually superceded by the analytic result, given the huge 
number of Feynman diagrams to evaluate by tedious recurrence relations, having information on 
the final results from an independent calculation provided an important crosscheck on the full 
n-dependent expressions. This is particularly the case when computations of gauge invariant 
quantities are simplified by choosing to work in the Feynman gauge. This substantially reduces 
the large number of integration by parts and hence the computation time is significantly smaller 
than, say, in an arbitary linear covariant gauge. In other words working in a specific gauge 
means that the internal strong check of observing the gauge parameter cancellation for a gauge 
independent quantity is crucially absent. 

Given the potential determination of analytic anomalous dimensions for the transversity 
operator in the foreseeable future, it is the purpose of this article to repeat the approach of [14] 
and provide the MS anomalous dimensions for fixed moments at three loops up to and including 
moment n = 8. This builds on the low moment results of |171ll8j where the anomalous dimensions 
of the tensor current and second moment were determined. More recently the results for moments 
n = 3 and 4 were provided in |19j . In |18l [T9] the primary aim was to provide the finite parts 
of a specific Green's function in order to aid lattice measurements of the same quantity. In 
particular those measurements had to match onto the ultraviolet part of the Green's function 
and the provision of the answer at three loops was necessary to help make the extraction of lattice 
results as precise as possible. However, the work also required performing the renormalization in 
the lattice renormalization scheme, known as the modified regularization invariant (RI') scheme, 
[20\ [2T] . It has a continuum definition which is discussed at length in [22\ ITT] . As the anomalous 
dimension of a gauge invariant operator is gauge dependent in a mass dependent scheme, the 
computation had to be performed in an arbitrary linear covariant gauge. Though for practical 
reasons the lattice calculations were performed in the Landau gauge. Therefore, as a second 
thread to this article, we will also determine the transversity anomalous dimensions in the RI' 
scheme in an arbitrary linear covariant gauge. Though for computational reasons this will be 
restricted to n < 7. 

First, we discuss the basic properties of the transversity operator and outline our computa- 
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tional strategy. The operator is defined by, [HI [9| [TU]. 

0m...i^n ^ Si^a'^'^D"^ ...D^^iP (1) 

where a^'^ = ^['~f^,j'^] and is the covariant derivative involving the couphng constant g. The 
operation 5 denotes symmetrization of the Lorentz indices {i^i, . . . , Vn} as well as ensuring that 
the operator is traceless according to the rules 

^ Q (i > 2) 

r?.,.,C"^'^^-"-"^-"" =0. (2) 

To renormalize ([T]) we follow a procedure similar to [171 118j where the operator is inserted at 
zero momentum into a quark 2-point function, {ip{p)O{0)ip{—p)), where p is the momentum. 
In order to apply the Mincer algorithm, [23], written in the symbolic manipulation language 
Form, [241 [25] . one needs to saturate the Lorentz indices with the appropriate tensor. This 
is because the Mincer formalism can only be applied to massless three loop Lorentz scalar 
2-point Feynman integrals, [23|. Here, since we are not interested in the finite part of this 
Green's function, we merely multiply it by a Lorentz tensor which has the same symmetry and 
tracelessness properties as the original operator itself. This leads to an immediate algebraic 
simplification. When the fully symmetrized and traceless operator is inserted into the Green's 
function, there is a part involving products of the tensor 77^*^. These derive from ensuring the 
tensor is overall traceless. However, when these terms multiply the projection tensor such terms 
will give zero. Therefore, in constructing the Feynman rules for the operator insertion for the 
current calculation, one needs only to consider the part of the operator which is independent 
of the r]^^'^ tensors. In other words, the object ^pa'^^'^^ . . . D'^"^'ijj. For high moments, this 
represents a huge reduction in work such as the tedious but automatic derivation of the full 
Feynman rules, which can involve a significantly large number of terms. Instead the main work 
is in the construction of the projection tensor. However, this is achieved computer algebraically 
by writing down the complete set of independent objects built from one a^'^ tensor, together 
with the appropriate numbers of r]^'^ tensors and momenta p'^ such that the number of free 
Lorentz indices equates with that of the original operator. These independent tensors are then 
symmetrized automatically with respect to the indices {vi, . . . , and the arbitrary coefficients 
chosen so that the overall projection tensor is traceless according to ([2]). As was indicated in 
\18\ I19j . there are three independent projections and therefore for the determination of the 
anomalous dimensions one needs only to select one of these for the projection procedure. It is 
preferable to choose the most algebraically compact one to minimize computation time. We note 
that this procedure is in contrast to the fixed moment strategy of [14] where the null vector 
was introduced. However, as was discussed in |141l26j. one needs to handle the resulting integrals 
using a different projection technique. Finally, we note that our procedure automatically ensures 
that there is no operator mixing, [19j. The Feynman diagrams required for {'ip{p)O{0)'ip{—p)) are 
generated automatically with the Qgraf package, |27j . and converted to Form input notation 
prior to the application of the three loop Mincer algorithm. There are 3 one loop, 37 two loop 
and 684 three loop Feynman diagrams and throughout we have used dimensional regularization 
in d = 4 — 2e dimensions. 

We now record that the three loop MS transversity anomalous dimensions, 7^"^(o), for n = 5, 
6 and 7 are 

7(5) (a) = —Cpa + [189515C7A - AmACp - rdSlOTpNf] --|— 
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(190836000C(3) + 1975309075) 0% 
- (572508000C(3) + 325464235) CaCf 
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- (1192320000C(3) + 511395100) CATpNf 

+ (381672000C(3) - 254723696) C| 

+ (1192320000C(3) - 989903260) CpTpNf 



83718800T§^ Nf 



12150000 



(3) 



y C^a + [204770C^ - 42129C7^ - 88810rF7V/] 



+ 



(707616000C(3) + 7527909825) Cl 

- (2122848000C(3) + 1373507730) CaCf 

- (4626720000C(3) + 1841332000) CATpNf 
(1415232000C(3) - 684744816) 
(4626720000C(3) - 3910683210) CfTeNj 



+ 
+ 



320975800r|.iV/ 



42525000 



+ 0(a* 



(4) 



and 
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"35" 

+ 



Cpa + [75266555Ca - 15484767Cir - 33149830rpA(/-] 

(3517994592000C(3) + 38365845513450) C\ 

- (10553983776000C(3) + 5978407701105) CaCf 

- (24084527040000C(3) + 9039144860900) CaTfNj 
(7035989184000C(3) - 4192441946262) C| 
(24084527040000C(3) - 20698675427220) CpTrNf 



_CW_ 
2315250 



+ 
+ 



1651311191600r|iV/ 



204205050000 



+ 0(a^ 
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where a = (7^/(167r^), (^(n) is the Riemann zeta function, is the number of quark flavours 
and the colour group Casimirs are given by 



Tr ( 'Y'^T^ 



T^T" = Cpl 



jacd jhcd 



CaS 



ab 



(6) 



and are the generators for the colour group with structure constants We have also 

calculated the n = 8 moment but modified the method used to find ([3]), (|4|) and ([5]). In order to 
reduce the computation time we performed the calculation in the Feynman gauge. Whilst this 
significantly reduces the number of integration by parts needed to be carried out by the Mincer 
algorithm, for example, the check of gauge invariance of the final result in the MS scheme is 
removed. However, as a partial check on the automatic construction of the n = 8 Feynman 
rules and projection tensor, we have performed the two loop calculation in an arbitrary linear 
covariant gauge and reconstructed the gauge independent result of fSl O [ini [HI [IS]. Thus, at 
three loops in the MS scheme we have 



551 

+ 



Cpa + [1270588235CA - 251839827(7^ - 568470280rpiV/] 

(57548150352000C(3) + 651153115163775) C\ 
- (172644451056000C(3) + 110529966326535) CaCf 



Cpa^ 



37044000 
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- (411490679040000C(3) + 148380288276500) CaTfNj 

+ (115096300704000C(3) - 55777651074312) C], 

+ (411490679040000C(3) - 356756758487220) CpTpNf 



27902636165600T|iV/ 



3267280800000 



+ 0(a^ 



(7) 



As further checks on our results, ([3]), ([1]), ([5]) and ([7]), we note that for MS all the two 
loop expressions agree with those of [U O [ini ttH US], when they are evaluated for specific n. 
Moreover, ([3|), @ and ^ are clearly independent of the linear covariant gauge fixing parameter, 
a. Further, the coefficients of the leading order large Nf term at three loops, corresponding to 
the CpTj^Nf term, agrees with the analytic evaluation at 0{1/Nf) for arbitrary n given in [18]. 
A final check is provided by the fact that the triple and double poles in e of the resulting operator 
renormalization constant are reproduced precisely in agreement with the renormalization group 
equation prediction for all four cases. This is important since we followed the procedure of [28] 
for renormalizing operators in automatic symbolic manipulation programmes. For the specific 
case of the Lie group SU{3) we have 



7(^)(a) 



SU(3) 368 

= 1^° - 2 IlliS^lW;- 15389391 — 

- 188367300iV/ + 8942400000C(3)iV/ + 12843253410A[/- 



954180000C(3) - 144207743479 



82012500 



+ 0{a' 



(8) 



5(7(3) 136 a 

= —a - 2 [U405Nf - 558138] 



15 ' ' 10125 

240731850iV/ + 11566800000C(3)A^/ + 16107360420A?^ 

- 1179360000C(3) - 183119500163 



95681250 



+ O(a') (9) 



and 



5(7(3) 



SU{3) 



344 

a - 2 [16574915iV/ 

35 



205153309] 



3472875 

206413898950iV/ + 10035219600000C(3)Ar/ 
+ 13678917121415/Vy - 977220720000C(3) 



- 156962874344971 



76576893750 



+ O(a^) 



ii^a - [284235140Ar. - 3475978269] 
105 ^ ^ 27783000 

- 20926977124200iV/ + 1028726697600000C(3)iV/ 
+ 1381224814218690Ay - 95913583920000C(3) 



15957293707773841 



7351381800000 



+ 0{a^ 



(10) 



(11) 



where Cp = 4/3, Ca = 3 and Tp = 1/2. 
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Finally, we note that the results analogous to ([l]) and ^ in the RI' scheme are 



92^ 



+ 



(30825^2 + 92475a + 174069o) Ca - 166696C^ - 676560TpiV/ 



27000 

,2 



(l94197500a'^ + 1854279000a^ - 583200000C(3)a2 + 8993896875a^ 
- 6026400000C(3)q + 30074295375q - 37353312000C(3) 
+ 356401468700) 

+ (91239750a^ - 2099569500^ - 4997987400a + 1076976000C(3) 

- 60979980560) CaCf 
- (17262000000^ - 1555200000C(3)o + 10041363000o + 17858880000C(3) 

+ 253330505600) CATpiV/ 
+ (1289803200O + 27164160000C(3) - 22363266560) C^TpiV^ 



+ (10686816000C(3) - 7132263488) C|, + 41629683200r|A(^2 



340200000 



+ 0(a' 



(12) 



and 



—Cpa 
5 



(330750^ + 992250 + 195791o) Ca - 168516Ci7' - 769360rpA^/j 

(ll9070000o^ + 11600280000^ - 364500000C(3)o2 + 56219546250^ 
- 3766500000C(3)a + 18833983125a - 22639824000C(3) 
+ 231206558900) 

+ (494370000^ - 1764774000^ - 3356146800a - 1150848000C(3) 

- 40247461840) CaCf 
- (10584000000^ - 972000000C(3)o + 6258495000a + 11741760000C(3) 

+ 165502535600) CaTp TV/ 
+ (866102400a + 18040320000C(3) - 14470193920) CpTpNi 

+ (6469632000C(3) - 3130262016) C|, + 27562342400r|A^? .^..^^o^o 

J 194400000 

+ 0(a4) (13) 



258 
"35" 

+ 



+ 



(l2006225a2 + 36018675o + 73991771o) Ca 

C 

- 61939068CF - 292181680rpA^/] ^^^^ 

(15884235675000^^ + 15732651627000o^ - 4900921200000C(3) 

+ 761656331913750^ - 50642852400000C(3)o + 255327460003875a 
- 301268627856000C(3) + 3230428196219700) C\ 



o 
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+ (598161044250a^ - 2827913095350a^ - 47764761023700a 

- 30807568512000C(3) - 558295101833000) CaCf 
- (l4119320600000a2 _ 13069123200000C(3)a + 84542124471000a 

+ 164110847040000C(3) + 2326270368265200) CATpNf 
+ (12326389941600a + 254163329280000C(3) - 198131779069280) CfTfNj 
+ (84431870208000C(3) - 50309303355144) C|, 



+ 391057633881600riiV? + Oia^) . (14) 

^ ^ J 2450460600000 ^ ' ^ ' 

To summarize, this renormalization scheme is defined by ensuring that the finite part of the 
renormahzed Green's function {jp(p)O{0)tlj{—p)) , multiphed by the projector, is given purely by 
its tree value only, [2D1 EH [HI [ISl [TO] . Further, it is important to note that in i^, ([H]) and (HID 
the variables a and a are to be regarded as RI' quantities. The relation to the MS variables are 
given to three loops in [T7] for an arbitrary linear covariant gauge. Since our main motivation 
is to provide the MS anomalous dimensions, the n = 8 moment in the RI' scheme is clearly not 
available since the corresponding MS computation was restricted to the Feynman gauge. 

We conclude by noting that we have provided higher moments of the transversity operator at 
three loops in both the MS and RI' schemes. Together with the earlier results of 120 } 1211 [T8l ll9j. 
there are now eight fixed moment MS anomalous dimensions available prior to an explicit n- 
dependent computation. This is similar to the situation with regard to the flavour non-singlet 
twist-2 operator where there were seven fixed moment anomalous dimensions available prior 
to the provision of the full n-dependent expression. Whilst it is still in principle possible to 
compute even higher moments using the method we have discussed here, we believe we have 
reached a computational limit beyond which it is not viable to proceed. For instance, for the 
three loop re = 8 Feynman gauge calculation, the necessary Feynman rules took around 36 hours 
to be generated electronically on a dual opteron 64 bit SMP machine (2 GHz), resulting in just 
under 5.5 x 10^ terms. The former is an order of magnitude larger than the time required for 
the re = 7 case. 
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